Bloch oscillations of Path-Entangled Photons 
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We show that when photons in iV-particle path entangled |iV, 0) + |0, N) state undergo Bloch 
oscillations, they exhibit a periodic transition between spatially bunched and antibunched states. 
The transition occurs even when the photons are well separated in space. We study the scaling of 
the bunching-antibunching period, and show it is proportional to 1/N. 



When electrons in crystalline potentials are subjected 
to uniform external fields, classical mechanics predicts 
that they will exhibit Olmiic transport. Remarkably, 
in 1929 Bloch predicted that the quantum coherence 
properties of the electrons prevent their transport P, Q ■ 
He showed that the electrons dynamically localize and 
undergo periodic oscillations in space. Bloch oscilla- 
tions (BOs) manifest the wave properties of the elec- 
trons, and therefore appear in other systems of waves 
in tilted periodic potentials. BOs were observed for elec- 
tronic wavepackets in semiconductor supperlattices Q, 
matter waves in optical lattices Q and light waves in 
tilted waveguide lattices @|and in periodic dielectric 
systems [7J- 

In optics, BOs manifest the classical wave properties 
of light, and not its quantum (particle) nature. Recently, 
quantum properties of light propagating in periodic lat- 
tices of identical waveguides have been studied, predict- 
ing the emergence of nontrivial photon correlations 0, [j| . 
Nonclassical correlations between photon pairs were ex- 
perimentally observed in periodic lattices [lOj, while the 
effect of disorder was studied in jllj. BOs of a single pho- 
ton in tilted lattices were shown to follow the dynamics 
of coherent states [12j. Nonclassical features of BOs of 
photons in a two-band model were studied by Longhi, 
who showed that the probability to detect photon pairs 
in different bands oscillates nonclassically |l3l |. 

In this paper we study theoretically the propagation 
of spatially entangled states in waveguide lattices which 
exhibit Bloch oscillations. We consider light fields ini- 
tiated in a superposition of N photons in site // or in 

site |V) - ^ (|A%, |0)„, +e-^ |0> M , \N) u ,). Such 
superpositions, coined NOON states, exhibit fascinating 
quantum interference properties. NOON states are con- 
sidered the optimal quantum states of light for quantum 
meteorology applications such as quantum lithography 
and quantum imaging Here we show that when 

NOON states undergo BOs, the nature of the correlations 
between the photons oscillate between spatially bunched 
and antibunched states. We find that the period of the 
oscillations is inversely proportional to the photon num- 
ber N, resembling the X/N oscillations of NOON states 
in Mach-Zehnder interferometers. Interestingly, the os- 
cillation period is also inversely proportional to the initial 
separation of the two input sites — v' . A unique feature 



of the NOON state BOs is that the transition between the 
bunched and antibunched states can happen even when 
the photons are separated by many lattice sites. 

We consider the simplest waveguide structure which 
exhibits BOs, a one-dimensional lattice of single mode 
waveguides which are evanescently coupled. The prop- 
agation in the lattice is determined by two parameters: 
the phase accumulation rate in the waveguides (the prop- 
agation constant) and the tunneling rate between neigh- 
boring sites (the coupling constant) [15]. The propaga- 
tion of the fields in waveguide lattices is described by the 
tight-binding model, and was used to demonstrate many 
optical analogues of solid-state phenomena [161 [l7j ■ BOs 
are observed when the coupling constants between all 
the waveguides are identical and the propagation con- 
stants depend linearly on the waveguide position Q . 
To study the propagation of nonclassical light in such 
a structures we quantize the fields in the lattice. Since 
each of the waveguides supports a single mode, the field 
in waveguide [i is represented by the bosonic creation and 
annihilation operators a'^ and a M , which satisfy the com- 



mutation relations 



The operators evolve 



according to the Heisenberg equations 



i— - = uBa 

oz 



C 



(1) 



Here z is the spatial coordinate along the propagation 
axis, C is the coupling constant and B is the difference 
in the propagation constants of neighboring sites. The 
evolution of the creation and annihilation operators is 
calculated using the Green function U^^^z) of Eq. ([1}, 
a ji( z ) = Eu' U fl , fl '(z)a^,(z = 0) The unitary trans- 
formation f/ju,^' (z) describes the amplitude for the tran- 
sition of a single photon from waveguide fi to waveguide 
fj,' . The Green function of Eq. (QJ is given by 0, [l8| : 

EW(*) = eW-rieiW+rij^ (<§sin(Bz/2)) , 

(2) 

where J M (x) is the /ith Bessel function of the first kind. 
Since any input state can be expressed with the creation 
operators and the vacuum state |0), the evolution of 
nonclassical states along the lattice can be calculated us- 
ing Eq. ([2]) . The probability to locate at site /i a photon 
that is injected into the lattice at site // = is given 
by the photon density — (a^a^) = \U t 
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Figure 1: (Color online) (a) The photon density = 
(aLa M ) for a single photon initiated at the waveguide y! = 0. 
The photon is mostly localized in two narrow peaks at the 
edges of the distribution. The path of each peak has a period 
As, and marks a branch of the Bloch oscillation (red dashed 
line marks the right branch), (b) The photon density n M (z) 
for a NOON state input with N = 1 coupled to waveguides 
fj! — and v' = 1. The photon is mainly located at the right 
branch of the oscillation, (c) The photon density n^(z) for a 
NOON input state with N = 2 coupled to waveguides y! — 
and i/' = 1 The photons from y! = and i/ = 1 inputs add 
up incoherently, showing double-branch oscillations. 



is depicted in Fig. 1(a). The photon exhibits BOs: it 
spreads across the lattice by coupling from one waveg- 
uide to its neighbors in a pattern characterized by two 
peaks at the two edges of the distribution. Each peak 
covers approximately four waveguides, and and its loca- 
tion oscillates along the propagation axis with a period 
Xb = 2tt/B. The right and left paths of the two peaks 
mark the two branches of the BO. We note that such a 
double-branch pattern is not a special feature of single 
photons. Any state of light which is coupled to a sin- 
gle waveguide in the lattice will exhibit exactly the same 
photon density. However, when the light is coupled to 
more than one waveguide, the propagation of the pho- 
tons becomes state-dependent. Rai et al. have shown 
that a single photon which is coupled into the lattice in 
a superposition of several waveguides exhibits BOs like a 
coherent state [12]. Figure 1(b) shows the photon den- 
sity for a superposition of a single photon initiated in a 
superposition of two neighboring waveguides, with a rel- 
ative phase ip — 0. The two paths the photon can take, 
starting either from waveguide fi' = or from waveg- 
uide v' = 1, contribute coherently to the photon density 
n M = ^|?7 Al . Al ' = o + U^. U ' = i\ 2 . Because of this interference 
the photon oscillates in a single branch, exactly like a 
coherent beam. In contrast, when a NOON state with 
N > 1 is coupled to the lattice, the photon density is 
identical to the photon density obtained by two incoher- 
ent beams = f ll^^oO)! 2 + f \U fJ , y= i(z)\ 2 (Fig. 
1(c)). 

Nonclassical features of light are probed by correla- 
tions between the photons. We focus on the probability 
to detect p photons in waveguide /i and q — N — p pho- 
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Figure 2: (Color online) Bloch oscillations of NOON states 
with N — 2 coupled to two adjacent waveguides \ip) = 
-7j(|2) \0} 1 +e' l ' fi |0) |2) 1 ). (a) The multiple detection prob- 
ability T^Jv at several propagation distances, for ip — 0. At 
the beginning of the propagation the two photons exhibit an- 
tibunching and are located at the two different branches of 
the oscillations. As the photons approach the turning point 
(z — Ab/2), they bunch and are found with the highest prob- 
ability in the same branch, (b) Same as (a) for (p = ^. 
The photons show bunching-antibunching cycle, but in this 
case start the oscillation partially bunched, (c) Same as (a) 
and (b), for if — n. Here the photons start the bunching- 
antibunching cycle bunched, (d) The normalized coincidence 
rate 7' 1,1 '(z) as a function of the lattice length for <p — (blue 
solid line), ip = -| (green dash-dotted line), and <p = n (red 
dottde line). The coincidence rate is calculated between the 
positions of the central waveguide in each branch, showing 
oscillations with a period As. 



tons in waveguide v, T^fi 
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a NOON state coupled to waveguides //and v' , the mul- 
tiple detection probability is: 



r(P'i) — I 

L 2p\{N-p)\ 



(3) 



+ e l9(2) J,'- M (0 P ^-,(C) 9 



Where ( — AC /B sin (Bz/2), and 9{z) is given by: 



0(z) = <p+-(n + Bz)(i/-ii')N 



(4) 



Eq. ([3]) shows that two terms contribute to the multi- 
ple detection probability: the photons arrive either from 
the input waveguide // or from waveguide v' . Since the 
photons are indistinguishable, these two paths interfere. 
The phase between the two paths is proportional to Nz, 
indicating that the oscillation period scales like 1 /N (see 
below). 

In Fig. 2 we depict Tj^'v , the probability to detect one 
photon at waveguide \x and another photon at waveguide 
v, for a NOON state with N = 2. The left column shows 
the evolution of the probability for a NOON state with 
a phase ip = 0. At the beginning of the propagation 
(Bz <C 7r). the photons follow the same path as in a 
periodic array of identical waveguides 0, [ljj. At this 
stage the off-diagonal terms of the probability matrix 
r^^are much stronger than the diagonal terms, indi- 
cating that the photons exhibit antibunching: each pho- 
ton takes a different branch of the oscillation. However, 
during the expansion period of the BO, as the photons 
approach the turning point, the symmetry of the two- 
photon probability matrix changes significantly. The di- 
agonal terms of the matrix become more pronounced, 
i.e. there is a higher probability to find the two photons 
in the same branch. At the turning point z = ^f-, the 
photons bunch: the off-diagonal terms of the probabil- 
ity matrix vanish, indicating that the photons are never 
found simultaneously at the two different branches. Re- 
markably, even though the photons start the propaga- 
tion in spatially separated branches, at the turning point 
they bunch to one of the branches, with equal proba- 
bility. Beyond this point the photon density contracts 
back towards the origin waveguides. During this contrac- 
tion the pairs again switch to an antibunched state. We 
note that the bunching-antibunching transition happens 
when the two branches of the BO are spatially separated, 
whereas the bunching-antibunching transition predicted 
in binary lattices occurs only when the photons are in 
the same waveguide The cycle in the symmetry of 
the probability matrix is observed for any initial phase of 
the NOON state phase ip, as demonstrated Fig. 2(b) and 
2(c). The phase cp sets how bunched or antibunched the 
photons are at the beginning of the propagation, but the 
period of the cycle is phase- independent [see Eq. (QJ]. 
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Figure 3: (Color online) Bloch oscillations of NOON states 
with sub- As correlation-oscillation periods, (a) The multi- 
ple detection probability I^i, at several propagation dis- 
tances for the input state \i/>) = ^(|2)_ 1 |0) 1 + |0)_ 1 |2) 1 ). 
The photons exhibit bunching-antibunching oscillations (see 
text) with a period As/2. (b),(c) The multiple detection 



probability V 



{N/2, N/2) 



for a NOON state with N 



6 (b) 



and N = 10 (c), injected to adjacent waveguides \tp) = 
^(|A0ol°)i + PloWi)' The oscillations of the correla- 
tion matrix are much faster, hence the probability matrix is 
calculated for five lattice lengths close to the turning point 
z = As/2, (d) The normalized coincidence rate ^ N ^ 2 ' N ^ 2 \z) 
as a function of the lattice length for the above three cases. 
The period of the oscillations are As/2 [(a), blue solid line] 
As/3 [(b), green dashed-dotted] and As/5 [(c), red dotted 
line] . 
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The bunching-antibunching cycle described above can 
be realized experimentally by measuring the correlations 
at the output of lattices with identical parameters but 
with different propagation lengths. For each propagation 
length, the waveguide at the center of each oscillating 
branch (henceforth waveguides x and y) can be imaged 
on two photon-number resolving detectors. The proba- 
bility to detect p photons at waveguide x and q = N — p 
photons at waveguide y is proportional T^y . When a 
delay is introduced between the photons that are injected 
to waveguide // and the photons injected to waveguide 
u', the photons become distinguishable, as in the Hong- 
Ou-Mandel (HOM) experiment [2fj. This corresponds to 
replacing the NOON state with a mixed state of N pho- 
tons in either one of the two input waveguides. The ratio 
of the detection probabilities for the NOON and mixed 
states is given by 



Jp<q) 



-n(p.g) 

1 x,y 



1 jVj 

2 p\(N-p)\ 



I J p T q I 2 

\ J v'-x J u<-y\ 



(5) 

Figure 2(d) shows ^ 1 > 1 > as a function of the lattice 
length, for NOON states with N = 2 and <p = 0,tt/2, tt. 
When 'y( 1 > 1 J = 0, the photons are bunched and are never 
found in the two different branches of the BO; scanning 
the delay between the input ports of the lattice will yield 
a HOM dip. When 'y( 1 > 1 ) = 2, the photons are anti- 
bunched, and a delay scan will result in a HOM peak [2l| . 
Figure 2(d) clearly shows that the bunching-antibunching 
oscillations have a period of As. 

We next study input states which exhibit correlation 
oscillations with shorter periods. Eq. ^ suggests that 
the period of the oscillations in the correlation properties 
depends on the spacing between the input waveguides 
and on the number of photons in the NOON state. Fig- 
ure 3 shows several examples of correlation oscillations 
with periods shorter than Xb- In Fig. 3a we show the 
propagation for a NOON state with N = 2, where the in- 
put sites are separated by one waveguide. In this case the 
photons exhibit a bunching-antibunching transition with 
a different spatial symmetry [^| . The correlation map os- 
cillates between a state in which the peaks are highest at 
the corners of the correlations matrix, to a case in which 
the highest probability is between the corners. The os- 
cillation period is Ab/2, twice the period observed for 
a NOON state input with adjacent waveguides. Finally 
we calculate r^ /2,JV/2) for NOON states with N = 6 
[Fig. 3(b)] and N = 10 [Fig. 3(c)], with adjacent input 
waveguides. The oscillation period is indeed 2Xb/N, as 
predicted by Eq. (£[]). Within one oscillation of the single 
photon density, the TV-photon distribution switches N/2 
times from all the photons in the same branch to photons 
divided equally between the two branches. 

In conclusion, we studied the propagation of photonic 



NOON states in waveguide lattices which exhibit Bloch 
oscillations. We found that while the photon density os- 
cillates in the Bloch frequency, the multiple detection 
probability oscillates at higher frequencies. These oscil- 
lations indicate that the photons show a transition from a 
bunched to antibunched states, with a period that scales 
as 1/N. By carefully designing the parameters of the 
Bloch lattice this oscillatory transition can be used to 
distribute bunched and antibunched states of light in 
an integrated and thus robust manner. To experimen- 
tally observe the bunching-antibunching transition, we 
propose to perform a Hong-Ou-Mandel measurement be- 
tween two waveguides at the two branches of oscillations 
using photon-number resolving detectors. We predict os- 
cillations between a HOM dip and peak as a function of 
the propagation distance in the lattice. Recent progress 
in waveguide lattice fabrication [22I |23| . photon num- 
ber resolving detectors and photonic NOON state sources 
[2^ | . make such measurements in reach. 
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